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P. Holmes.  Research  Activities  January  1.  1988  -  September  31.  1988 

1.  Dynamics  of  0(2)  Symmetric  Systems 

J.  Guckenheimer .  D.  Armbruster  (postdoc)  and  S.  Campbell  (grad  student)  and 
I  have  continued  our  work  on  the  global  dynamics  and  bifurcations  of  0(2) 
symmetric  ODEs.  Such  systems  are  obtained  as  finite  dimensional  projections  or 
reductions  of  spatially  translation-  and  reflection-invariant  PDE’s,  for 
example.  In  1987/88,  partially  supported  by  this  grant,  we  provided  a  complete 
analysis  of  heteroclinic  cycles  and  modulated  travelling  waves  in  two  mode 
(k:2k)  interacting  systems  [1,2].  In  particular,  we  pointed  out  that 
heteroclinic  cycles  are  structurally  stable  features  in  such  systems  [cf.3]. 

The  original  suggestion  that  1nte**estlng  global  d3rna.'T>ical  beli^tvior  arises 
in  0(2)  symmetric  problems  came  from  work  on  boundary  layers  [4],  where  the 
heteroclinic  cycles  correspond  to  bursting  behavior.  In  that  system  five  or 
more  "modes"  are  involved  and  we  are  therefore  now  studying  a  three  mode 
(k:2k:4k)  interaction,  which  appears  to  capture  the  essence  of  a  crucial 
interaction  mechanism  in  the  five  mode  boundary  layer  model.  The  system  of 
ODEs  in  six  dimensional,  having  3  complex  (Fourier)  modes,  and  so  presents  a 
formidable  analytical  challenge.  Armbruster  and  Campbell  have  performed 
extensive  numerical  simulations  and  uncovered  pervasive  families  of  modulated 
travelling  waves.  We  are  currently  studying  various  types  of  limiting  systems, 
for  small  parameter  regimes  near  the  multiple  bifurcation  point  at  which  all 
three  modes  are  neutrally  stable,  as  well  as  attempting  to  approximate  Poincare 
return  maps  near  the  heteroclinic  cycles.  We  hope  to  continue  this  work  under 
a  new  grant. 

2.  Spatially  Complex  Equilibria  of  Buckled  Rods 

The  work  done  by  A.  Mielke  and  me  in  1987  has  been  published  [5]  and  these 
studies  are  now  complete.  Using  "Melnikov”  perturbation  methods  and  exploiting 
the  noncainonical  Hamiltonian  structure  of  the  equilibrium  equations  for  a 
hyperelastic  rod,  we  were  able  to  provide  explicit  computations  demonstrating 
the  existence  of  spatially  "chaotic"  equilibrium  shapes.  Such  solutions  are 
the  spatial  analogue  of  the  chaotic  dynamics  proven  to  occur  in  the  dynamical 
problems  of  a  rotating  rigid  body  subject  to  gravity  (cf.[6],  Kirchoffs’ 
kinetic  analogy). 
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1.  Averaiging  using  Jacobian  Elliptic  functions 

We  have,  together  with  graduate  student  Vincent  Coppola,  investigated  the 
dynamics  of  the  system 

(1)  x’ ’  +  a(T)  X  +  /3(t)  x^  =  fc  g{x,x’,t) 

in  which  t  =  et  is  slow  time.  The  idea  here  is  that  when  e  =  0  the  system 
reduces  to 

(2)  x’  ■  +  a(0)  X  +  j3(0)  x^  =  0 

which  has  a  solution  in  terms  of  Jacobian  elliptic  functions,  e.g. 

(3)  X  =  A  cn(<p.k) 

We  have  applied  the  method  of  averaging  to  first  and  second  orders  in  e 
to  the  system  (1).  This  involves  perturbing  off  of  Jacobian  elliptic 
functions,  rather  than  off  of  trignometric  functions  as  is  usually  done.  The 
resulting  equations  involve  integrals  of  elliptic  functions  which  are  evaluated 
using  a  program  written  in  the  computer  algebra  system  MACSYMA.  The  results 
have  been  applied  to  the  problem  of  finding  limit  cycles  in  the  autonomous 
version  of  eq.(l)  [1],  as  well  as  to  the  chaotic  behavior  of  the  eq. : 

3 

(■i)  x’ '  +  X  cos  T  +  X  =0 


The  origin  in  this  equation  is  saddle-like  for  half  of  its  period,  and  is 
center-like  for  the  other  half.  By  analyzing  the  slow  flow  associated  with 
eq. (4), 


k  1  w  « 
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we  have  been  able  to  prove  the  existence  of  chaotic  behavior  in  the  eq.(4). 

2.  Computer  Algebra  Treatment  of  the  Elliptic  Restricted  3  Body  Porblem 

We  have,  together  with  graduate  student  Bhaskar  Viswanadham,  investigated 
the  stability  of  the  equilibrium  solution  in  the  elliptic  restricted  problem 

of  three  bodies.  This  problem  is  governed  by  the  following  d.e.’s- 


(2.1) 


x”  -  2  y*  -  g  hg  X  =  0 


(2.2)  y’*  +  2  X-  -  g  hj  y  =  0 

where  g  =  1/(1  +  e  cos  f)  and 

^.2  =  I  i  "  3p(l-p)]^/2j 


in  which  primes  represent  differentiation  with  respect  to  the  true  anomaly  f 
(the  independent  variable).  Eqs.(2)  Involve  two  parameters  p  (mass  ratio  of 
the  smaller  primary  to  the  sum  of  the  primaries)  and  e  (eccentricity  of  the 
orbit  of  the  primaries).  For  given  p.e,  the  equilibrium  is  stable  if  all 

solutions  to  (1)  are  bounded,  and  unstable  if  an  unbounded  solution  exists. 

Although  numerous  papers  have  been  written  since  the  mid-1960’s  on  this 
problem,  we  have  gone  far  further  in  the  perturbation  analysis  than  previous 
investigators  by  the  use  computer  algebra  and  an  efficient  perturbation 
algorithm.  We  obtained  45  terms  in  the  perturbation  series  for  the  transition 
curves  separating  regions  of  stability  in  the  p-e  plane  from  regions  of 
instability,  with  the  result  that  studies  of  the  convergence  of  the  series  have 
been  made  possible.  It  was  fotind  that  one  of  the  series  exhibits  unusual 
behavior  in  that  it  can  be  shown  to  converge,  but  not  to  the  correct  (i.e., 
numerically-obtained)  solution  [2]. 
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F.Moon.  Research  Activities  January  1.  1988  -  September  31.  1988 

1.  Chaotic  Vibrations  of  a  Three-Well  Potential  Oscillator 

Together  with  postdoctoral  research  associate  G.  Li,  we  have  studied  the 
forced  vibrations  of  a  nonlinear  elastic  beam  having  three  stable  and  two 
unstable  equilibrium  positions.  The  responses  of  the  system  exhibit  chaotic 
behavior  in  certain  regions  in  the  driving  frequency  and  amplitude  plane.  With 
the  use  of  the  Melnikov  technique,  the  homoclinic  and  heterocliric  bifurcations 
of  the  system  were  calculated.  Numerical  simulations  of  the  stable  and 
unstable  manifolds.  Poincare  maps,  Lyapunov  exponents,  and  fractal  basin 
boundaries  were  performed.  The  results  indicate  that  both  homoclinic  and 
heteroclinic  bifurcations  are  necessary  conditions  for  chaos,  and  that  the 
final  steady  state  of  motion  for  many  initial  conditions  becomes  unpredictable 
if  at  least  one  homoclinic  or  heteroclinic  bifurcation  occurs.  An  experimental 
criterion  and  a  heuristic  formula  for  chaotic  motion  were  obtained  eind  compared 
with  the  bifurcation  criteria. 

2.  Chaotic  Nonplanar  Motions  of  a  Forced  Linearly  Elastic  Rod. 

Together  with  graduate  student  J.Cusumano  we  have  combined  a  geometrically 
exact  rod  theory  with  new  analytical  and  experimental  techniques  from  dynamical 
systems  theory  to  explain  complex  dynamical  phenomena  observed  in  a  thin 
prismatic  steel  rod  with  clamped-free  boundary  conditions.  The  support  of  the 
rod  is  harmonically  displaced  with  the  axis  of  the  displacement  aligned  with 
the  lateral  axis  of  symmetry  of  the  rod  so  that  one  would  expect  motions  to 
remain  planar.  However,  experiments  show  that  planar  motions  become  unstable 
when  the  forcing  frequency  is  close  to  a  resonant  frequency  of  the  system.  The 
resulting  motions  are  nonplanar  and  chaotic. 

A  previously  undiscovered  family  of  asymmetric  bending-torsion  nonlinear 
modes  have  been  found  experimentally.  Other  phenomena  discovered  in  the  system 
include  dynamic  two-well  behavior  and  energy  cascading  from  high  frequency 
input  to  low  frequency  response.  The  fractal  dimension  of  the  attracting  sets 
have  been  estimated  directly  from  experimentally  obtained  scalar  time  series. 

The  resulting  dimension  estimates  are  below  5,  implying  that  it  may  be 
possible  to  model  the  dynamics  of  the  rod  with  as  few  as  two  degrees  of 
freedom.  A  geometrically  exact,  linear  elastic  rod  theory  is  developed,  and 
used  to  derive  a  model  system  of  two  nonlinear ly  coupled  pxie’s,  which  is  then 
further  approximated  by  two  ode's  by  the  assumed-modes  method.  A  family  of 
nonlinear  modes  analogous  to  those  found  experimentally  in  the  rod  were 
discovered  numerically.  The  behavior  of  the  model  has  been  favorably  compared 
to  the  experimental  data. 
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